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Abstract 

The correlation functions of the spin-^ XXZ spin chain in the ground state are expressed in the 
form of the multiple integrals. For —1 < A < 1, they were obtained by Jimbo and Miwa in 1996. 
Especially the next nearest-neighbour correlation functions are given as certain three-dimensional 
integrals. We shall show these integrals can be reduced to one-dimensional ones and thereby 
evaluate the values of the next nearest-neighbor correlation functions. We have also found that 
the remaining one-dimensinal integrals can be evaluated analytically, when v = cos _1 (A)/-7r is a 
rational number. 
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The spin- 1/2 XX Z chain is one of the fundamental models in the study of the low- 
dimensional magnetism. The Hamiltonian is given by 



(i) 



j=-oo 



where S = cr/2 and o are Pauli matrices. The model can be solved by Bethe ansatz method 
and diversephysical properties have been investigated with varying the anisotropy parameter 
A |ll uj y|, |4| . Especially, in the region — 1 < A < 1, the ground state is critically disordered 
and the excitation spectrum is gapless. Then the long-distance asymptotics of two-point the 
correlation functions such as {S°/S°/ +k )k-^\ <^ = x,y,z are shown to decay as a power law 
via field theoretical approach (see, for example, and the references therein). However, if 
possible, it is more desirable 

(1) to calculate (S^S?^) for finite k first and 

(2) to derive its asymptotic behavior exactly. 

Unfortunately such a program has not been achieved except for the A = case 

00. 

In 1996, Jimbo and Miwa [^] obtained the multiple integral representation for the arbi- 
trary correlation functions of the XX Z chain for — 1 < A < 1. For example, the Emptiness 
Formation Probability (EFP) {J defined by 



Pin) 



(2) 



has the integral representation 
P(n) = 
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where the parameter v is related to the anisotropy A as 
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22j. Particularly, in the isotropic limit 
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A — > l{y — > 0), the general method to evaluate the multiple integral was recently developed 
by Boos and Korepin 



It is of some importance to note that P(2) and P(3) are related to the nearest and next 
nearest-neighbor correlation functions : 



(S*S* +1 ) = P{2) - 1/4, (5) 
(S*S* +2 ) = 2(P(3) - P(2) + 1/8). (6) 

One of our purposes in this letter is to evaluate (SjSj +2 ) through the integral representations 
of P(3) and P{2). Similarly the nearest and the next nearest-neighbor transverse correlation 
functions have the integral representations, 

, x x , 1 f°° dxi f°° dx 2 sinh(x 2 — x{) sinh ((x% + 'm/2) v) sinh ((x 2 — i7r/2) v) 
\ b j b j+i) - ~ — 



2v 2tx J_ 00 2tt sinh ((x 2 - x-i) v) cosh 2 x\ cosh 2 x 2 
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and 



k=i j - 



1 -A- f°° dx k sinh 3 ^ ((x k + m/2) v) sinh^ 1 ((x k - m/2) v) 



2ix cosh 3 Xk 



^ sinh(x2 — xi) sinh(a;3 — x\) sinh(x3 — x 2 ) ^\ 

sinh {{x 2 — X\) v) sinh ((£3 — x\) v) sinh ((0:3 — x 2 — lit) v) ' 

respectively. 

It was already shown by Jimbo and Miwa that the two-dimensional integral for P(2) 
reduces to the one- dimensional one for arbitrary v as 

„, . 1 1 d f f°° sinhfl - u)w , 1 

P(2) = - + — ^ — {sin™/ — -± r^-dtc , 9 

2 27r^ sin 7r^ az/ [ > /_ 00 sinh u; cosh z/w J 

which leads to the one-dimensional integral representation of (SjSj +1 ) via the relation 
The result coincides with a different derivation of (SjSj +1 ) from the ground-state energy 
per site e : 

de 1 <9e 



where 



dA 7rsin7rz/ <9z/ 

l^cot7rz/ r°° dw sinh(l — v)w 1 r°° dw wcoshw 
4 2ir J_ 00 smhw coshz/w 27r 2 7_ 00 sinh w (cosh z/w;) 2 ' 



A sm7ri> f smh(l — v)w . , 

e = — / — dw. (11) 
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Similarly we can get the one-dimensional integral representation for (SjSj +1 ) as 
1 



1 



oo 



dw; sinh(l — v)w cos tiv 



du> w cosh io 



(12) 



47rsin7rz/ J_ OQ sinhw coshz/w 47T 2 sinh w (cosh z/w) 2 

Thus, for the nearest-neighbor correlation functions, we have known the two-dimensional 
integrals by Jimbo-Miwa formula can be reduced to one-dimensional ones. The main result 
of this letter is that three-dimensional integrals for P(3) and (SjSj +2 ) can also be reduced 
to one- dimensional integrals. In other words, we have succeeded in performing the integrals 
for P(3) and (SjSj +2 ) twice. Our results are 

! 1 — cos27rz/ d f cosh3z/x 
dv \ (sinh , 
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Here 5 (0 < |5| < tt) should take some non-zero real value, which is introduced to avoid the 
singularity at the origin. We, however, note that the singular term may be subtracted from 
the integrand in principle, as its residue vanishes. Or alternatively, by applying the Fourier 
transform, we can express the one-dimensional representations (|13|) and (|14j) as 



i>(3) = - 2 + 
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Note that in the representations (|15p and (j!6|) . there are no singularities at the origin. This 
is in the similar situation as (fTU|) and (fT2"|) . Combining © and (fT5j) . we can also write down 
the one-dimensional integral representation for (SjSj +2 ) through the relation (jUJ), 



(SjSj +2 ) 



+ 



dw sinh(l — v)w 
sinh w cosh z/w 
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Now let us discuss some properties of the obtained one- dimensional representations |T6|) 
and (fTTj) . At a first glance, some integrands in (fTT)j) and (fTTj) are divergent at z/ = 1/2 and also 
in the limit v — > and z/ — ► 1, due to the the factor 1/ sin27rzA However, by investigating 
the integrands carefully, we have found these singular terms cancel each other, therefore 
yielding definite finite values for the correlation functions : 



1/2: 



It" 



v -> : 



/ QX QX \ I QZ QZ \ 



14 3 

-In2 + 4C(3), 

12 3 4 SV ; ' 



(18) 
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1/ -> 1 : 
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(20) 



It is especially intriguing to observe that the v — > limit ()19|1 reproduces the known result 
by one of the authors in 2^|. 

More generally, we have found when v takes a rational value, the one-dimensional integrals 
can be evaluated analytically. Some of our explicit results are summarized in Table 1. One 
can see that the correlation functions are in general a polynomial of 1/tt. Particularly, when 
v = 1/3, we expect all the correlation functions are given solely as a single rational number 
(c.f. I13L ll^). In Fig. and Fig. |2l plotted are the numerical values of the nearest and 
the next nearest-neighbor correlation functions calculated from the one-dimensional integral 
representations. For comparison, the analytical values in Table 1 are represented by the 
filled circles. 
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We have obtained (|T3j) and by generalizing the method developed by Boos and 
Korepin tl7 |. which allows us to calculate the multiple integrals for the correlation 
functions of the XXX model (A = 1). Below we outline the derivation of (|13|) and (|14|) 
quite briefly. The details of the calculations will be published in a separate paper jsjj. First 
we introduce the following convenient notation for P(3) : 

"°°- i/2 dAi 

-oo-i/2 



3 /-oo-i/2 j\ 

P ( 3 ) = IT / T-^(Ai, A 2 , A 3 )r 3 (A 1 , A 2 , A 3 ), (21) 



where 



and 



71-3 T7i<- 1..- sinh Ti (A,- — Afc) 
f/ 3 (A 1; A 2 , As) = - 3 ^f — ( 22 ) 



r 3 (Ai A 2 A 3 ) = ( qZl - ^ qZ2 - l ^ Z2 - ^ ~ ^ 2 (23) 
3 x ' 2 ' 3 8(z 2 - qzx)(z 3 - qz x ){z z - qz 2 ) 



with q = e 2niu ,Zi = e 2nuX % 
tion as Boos and Korepin 



i = 1,2,3). Then after similar but more complicated calcula- 



17j . we can transform the integrand T 3 (Ai, A 2 , A 3 ) into a certain 



canonical form without changing the value of the integral as 



where 



T C = P + — (24) 

Z2 - Zl 
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2q 8q 8q z Z\ 
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z 3 \ 8g 2 8g 8q 3 Zl r 

We find the first part of T c , namely Pq, can be integrated easily 



The second part P\/{z>i — Z\) can be integrated twice, namely with respect to A 3 and 
a = (A 2 + Ai)/2. Finally the remaining one- dimensional integration with respect to 
d = (A 2 — Ai)/2 together with (}2*Tj) gives the expression (fT3"]). 



Similarly for (SjSj +2 ), T 3 (Ai, A2, A3) in (|2T|) is replaced by 

r(\ \ \ \ - ( qZl ~ 1 ) 2 (g^2 ~ 1)(*2 ~ 1)(^3 - I) 2 , 9R s 

Ts(Al ' A2 ' " S q (z 2 -z 1 )(z 3 -z 1 )(z 3 -qz 2 ) ■ (28) 
As a corresponding canonical form, we have found 

T c = — ^— , (29) 



£2 - zi 



with 



(l + g)(3 + g 2 ) (3 - 2g + 3q 2 ) Zl (1 + qf 
n ^—5 1 „ H 



8g 2 8g 8g 3 ,2i 

3-q (l + q) Zl (l + gK-2 + g) 
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+ J L ^_ ii ^l_i_rji^ (30) 



4q 8q 8q 2 Zi 

l_ f l+_g _ (1 + _ 1+V 
23 I 2g 8 8g%! 

Again after integrating with respect to A 3 and a= (A 2 + Ai)/2, we arrive at the one- 
dimensional integral representation 1)140. 

In conclusion, we have shown the multiple integrals for the correlation functions of the 
XX Z chain at the critical region — 1 < A < 1, can be reduced to the one-dimensional ones in 
the case of the next nearest-neighbor correlation functions. This property will be generalized 
;o other higher-neighbor correlations as well as the correlations in the massive region (A > 1) 



2£,|2i|. In this respect, we like to refer to the recent work by Boos, Korepin and Smirnov 
ij , where they have shown the reducibility of the multiple integrals in the case of the XXX 
chain. 

In our future work, we are particularly interested in calculating the third-neighbor cor- 
relation functions (SjSj +3 ) and (SjSj +3 ) for gene ral A. For A = 1, they were recently 
calculated in [2^ from the multiple integrals by Boos-Korepin method . 
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the Scientific Research (B) No. 14340099 from the Ministry of Education, Culture, Sports, 
Science and Technology, Japan. GK and KS are supported by the JSPS research fellowships 
for young scientists. MS is supported by Grant-in-Aid for young scientists No. 14740228. 
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TABLE I: Some analytical values of the correlation functions when v takes rational values 
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